Collectively Enhanced Interactions in Solid-state Spin Qubits 
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We propose and analyze a technique to collectively enhance interactions between solid-state quan- 
tum registers composed from random networks of spin qubits. In such systems, disordered dipolar 
interactions generically result in localization. Here, we demonstrate the emergence of a single col- 
lective delocalized eigenmode as one turns on a transverse magnetic field. The interaction strength 
between this symmetric collective mode and a remote spin qubit is enhanced by square root of 
the number of spins participating in the delocalized mode. Mediated by such collective enhance- 
ment, long-range quantum logic between remote spin registers can occur at distances consistent 
with optical addressing. A specific implementation utilizing Nitrogen- Vacancy defects in diamond 
is discussed and the effects of decoherence are considered. 

PACS numbers: 03.67.-a, 76.30.Mi, 75.30.Hx 



Harnessing collective phenomena by utilizing ensem- 
bles of identical particles is a powerful tool, which has 
been exploited in effects ranging from superradiance to 
scattering suppression [l|. The coherent dynamics re- 
sulting from interactions with individual constituents of 
an ensemble are often too weak to be observed directly; 
however, as evidenced by experiments in systems such 
as Rydberg atoms [l-Q •> cavity QED [HL , atomic en- 
sembles @, IH, and solid state qubits [9[, collective en- 
hancement provides a natural route to overcoming this 
challenge. Here, we demonstrate that, for electronic spin 
quantum registers, such collective effects enable an ex- 
tended coherent coupling over large distances — an es- 
sential prerequisite for quantum information processing. 

Owing to favorable coherence properties, electronic 
spins associated with point-like defects in solid-state sys- 
tems have garnered significant recent interest as candi- 
dates for room-temperature quantum registers. Quan- 
tum control of such spins can be achieved using a com- 
bination of optical, magnetic and electric fields. While 
our considerations apply to a variety of electronic spin 
qubits here, we focus on the Nitrogen- Vacancy 

(NV) center in diamond. The NV center harbors an elec- 
tronic spin [S = 1), which can be optically initialized, 
coherently manipulated and read out on sub- wavelength 
scales [13l4l5] , These results have sparked several recent 
proposals which utilize networks of NV registers as the 
platform for a scalable quantum information processor 
[16|4l9|. However, for any spin qubit candidate, two cru- 
cial challenges remain to be addressed: 1) the weakness of 
the magnetic dipolar interactions on distances compati- 
ble with individual optical addressing and 2) the disorder 
in spin positioning due to inherent imperfections during 
defect creation. 



r = 20 nm 

FIG. 1: High-density NV spin ensemble distributed randomly 
within a sphere of diameter r, with an average distance a. The 
NV centers have three internal spin states that are split by a 
zero- field splitting and a Zeeman field. 



robust qubit with a collective, generally disordered spin- 
ensemble (Fig. [I]). If the spins behave in an aggregate 
fashion, such a qubit can produce a large state-dependent 
magnetic field, leading to enhanced long-range coupling 
between ensembles; this is reminiscent of tailored light- 
matter interactions achieved via atomic ensembles [20|. 
However, we note that quenched disorder naturally leads 
to localization in solid-state spin systems, implying that 
each eigenmode of the ensemble is composed of only a few 
spins. Here, we demonstrate the use of a uniform trans- 
verse magnetic field to overcome this issue. The applied 
field causes the symmetric I^-state 
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In this Letter, we present a novel approach to remote 
quantum logic which harnesses collectively enhanced in- 
teractions to overcome both of the above challenges. The 
key idea underlying our proposal is to associate a single, 



an approximate eigenstate of the Hamiltonian, thereby 
enabling us to harness it as a collective qubit. Moreover, 
we show that this particular state is largely insensitive 
to the underlying spin distribution and hence robust to 
effects of disorder. 

To be specific, we now describe our proposal in the 
context of NV diamond color centers. The largest energy 
scale in this system (A) is set by a combination of the 
zero field splitting (2.87 GHz) and a Zeeman field along 
the NV axis. We assume that this Zeeman field is suf- 
ficiently strong to ensure that the m s = — 1 spin state 
is sufficiently far detuned and hence does not contribute 
to the effective dynamics. Thus, the number of m s = 1 
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spins, m, is an approximately good quantum number and 
a perturbative description is justified. The second-largest 
energy scale arises due to the perturbation created by the 
transverse field £1. To gain a qualitative understanding, 
let us restrict ourselves to the analytically tractable case 
where m is either or 1. The effective Hamiltonian is, 
H r = -A|0X0| + y/NQ(\0){W\ + h.c), where the state 
|0) has all ensemble spins polarized into m s = 0, and 
the collective \W) state is fully symmetric with all spins 
sharing a single excitation, 



Hamiltonian including both on-site and interaction terms 
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Second order perturbation theory in y/Nfl/A yields 
H' r = -A|0X0| + J\W\W\, with J = NQ 2 /A. Including 
higher m manifolds merely leads to a renormalization of 
J, without changing this qualitative picture (so long as 
we are in the perturbative limit). This is equally true 
in the presence of dipolar interactions, provided that the 
energy scale J is larger than the characteristic strength 
of the dipolar interaction Vdd- Thus, even with these 
additional terms, the new eigenstates will still have sub- 
stantial overlap with the collective \W) state. This is in 
stark contrast to the situation without a transverse field, 
where strongly quenched disorder owing to random spin 
positions localizes all such eigenstates, even in three di- 
mensions. Furthermore, the dipolar interaction naturally 
ensures that collective states with different m values will 
have different energies, leading to a "blockade" -type sce- 
nario, where manifolds with m > 1 are energetically in- 
accessible 
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22| . This allows us to selectively drive the 
transition between |0) and \W) without populating any 
other collective states, provided that the external driving 
Qext is weaker than Vdd- This hierarchy of energy scales 
can be summarized as: A ^> J ^> Vdd ^> ^ext- 

Let us consider a three dimensional ensemble of TV = 
100 NV centers randomly distributed within a diame- 
ter r = 20 nm, as depicted in Fig. [U Such high density 
NV ensembles have been recently realized using long-time 
annealing of repeat-electron-irradiated diamond samples 
[2314261]. We will characterize our effective two- level sys- 
tem (m s = 0,1) using Pauli spin operators <j a . Being 
magnetic dipoles, NV centers interact with one another 
via long-range magnetic dipolar interactions (ignoring 
energy non-conserving terms which are suppressed by the 
NV center's zero field splitting), 



Vij = (l-3cos 2 ^ 
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where denotes the position, ji characterizes the mag- 
netic dipole moment, and 7?^ is the angle between the NV 
axis and the vector connecting sites and Yj. The total 



is then given by H = A/2 <y™ + ^E; + E*<j ^ • 
Now, let us turn to the enhanced coupling between 
an isolated NV defect (hereon termed "qubit") and the 
collective ensemble, separated by the distance R. We en- 
vision the ensemble to be initialized into the |0) state, 
while the NV qubit is initialized to the m s = 1 state. By 
ensuring that the qubit splitting is tuned resonant with 
only the \W) state, one finds that the effective dynam- 
ics are restricted to the single-excitation manifold of the 
combined qubit-ensemble system; to lowest order, these 
dynamics are governed by, 

#eff = >/Wc^(|l*,0X0 g , W\ + h.c), (3) 
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where N c characterizes the approximate number of spins 
participating in the \W) state and the notation |lg,0) 
refers to the combined state with the NV qubit being in 
m s = 1 and with the ensemble spins being in |0). Con- 
sistent with sub-wavelength techniques such as STED 
(R = 100 nm), we will assume that the NV qubit can be 
manipulated and read out independently of the ensemble 

To support the qualitative picture presented above, we 
now perform exact diagonalization of the full Hamilto- 
nian. In the majority of the numerics, we restrict our- 
selves to m < 2 excitations; however, we check the va- 
lidity of our results by including the m = 3 manifold for 
slightly smaller system sizes. For each eigenstate we 
calculate the collective enhancement factor, defined as 
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which essentially characterizes the number of ensemble 
spins participating in the eigenmode. As expected, in the 
absence of a transverse field, disorder localizes all eigen- 
states and as depicted in Fig. [2] (blue circles), N c <C N 
for all eigenstates. On the other hand, In the case of a 
moderate transverse field, Q « h x 100 MHz, one finds 
the existence of a single eigenstate with N c w 70 ~ N. 
While the specific details of this state depend on the mi- 
croscopic details (e.g., spin distribution within the en- 
semble and magnitude of the applied transverse field) , its 
collective nature is rather robust. In particular, as one 
varies the strength of the transverse field £7, there exists a 
large parameter regime where N c > 50 (Fig. [2]). The dips 
in N c are associated with resonance effects, which arise 
when other eigenstates become near-degenerate with the 
collective state. Finally, the decrease of N c for large val- 
ues of £1 signals the breakdown of perturbation theory as 
\fNVt/ A approaches unity. 

We now perform simulations of the combined qubit- 
ensemble system. As previously discussed, the system is 
initialized to |l g , 0) and the qubit splitting is tuned reso- 
nant with the energy of the collective mode; the resulting 
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FIG. 2: Comparison of the collective enhancement N c for 
Q = (blue) and Q = h x 110 MHz (yellow) within the single 
excitation manifold. In the latter case there is a collectively 
enhanced state with N c « 70, corresponding to an increase by 
more than one order of magnitude (A = /ix4 GHz). The inset 
shows the maximum value of N c depending on the transverse 
field strength Q. 



dynamics is evinced in Fig. [3] Interestingly, the probabil- 
ity of finding the qubit in the m s = 1 state, p q , exhibits 
collectively enhanced Rabi oscillations. The frequency 
of these oscillations is enhanced by nearly an order of 
magnitude relative to that expected for bare dipolar in- 
teractions between two individual NV qubits at a similar 
distance. The numerics also allow us to obtain the time 
required for an interaction-induced tt pulse, t n and from 
this, one can derive the effective distance R associated 
with i^eff • Surprisingly, in all cases, we observe that this 
distance corresponds not to R — r, but instead to the 
distance between the NV qubit and the center of the en- 
semble. We study the effects of putting the qubit closer 
to the ensemble by calculating the collectively enhanced 
coupling strength V c (as extracted from the numerically 
obtained t^). As shown in Fig. O we find that only for 
distances very close to the ensemble does the collective 
enhancement deviates from the asymptotic 1/R 3 scaling, 
e.g., the qubit is coupled to individual spins rather than 
to the entire ensemble. 

Experimental Realization and Decoherence. — Thus 
far, our discussion has assumed that both the NV qubit 
and the ensemble spins are perfectly decoupled from the 
environment. In any experimentally realistic scenario, 
however, there are two natural decoherence effects which 
will be present: spin dephasing and spin depolarization. 
We are particularly interested in the scaling of the de- 
coherence rates with TV, as this may adversely affect the 
scaling fidelity of our proposed long-range gates [Hj]. As 
the decoherence processes act locally on individual spins, 
we first calculate the decoherence rate for a single spin 
and multiply the result by N to obtain the rate for the 
collective state. For simplicity, we assume that the collec- 
tive state is the previously described \W) state in which 



a single excitation is shared among all N spins. 

First, let us consider the effects of spin dephasing. The 
worst-case scenario for such dephasing is given by the 
leaking out into non-symmetric states. Consequently, the 
error probability after a single T2 dephasing event on spin 
i is given by the probability to leave the \W) state, 
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where pt 2 is the single spin dephasing rate. For large TV, 
this result is essentially independent of N (after weigh- 
ing with the number of spins); therefore, the effect of T2 
processes on such a collective \W) state does not get en- 
hanced by system size and in fact, is only slightly worse 
than for a single spin, i.e., it can be expressed in terms 
of an effective dephasing time T| ff . 

Second, we consider the decoherence arising from 
phonon-induced spin depolarization processes (7\). 
Here, we must distinguish between processes which flip 
an ensemble spin from m s = 1 to m s = 0, and the re- 
verse. This asymmetry can easily be seen by noting that 
the \W) state has only one spin in m s = 1, while all 
other spins are in m s = 0. We denote the error probabil- 
ity associated with these two events as Pj^° and p^ 1 , 
respectively. For Pj^°, the state |0) with all ensemble 
spins in m s = is not affected at all, while the probabil- 
ity to flip from the \W) state into |0) is given by 



Pw^o=PT 1 ^o\(0\a i l ) \W}\ 
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which is again independent of the size of the ensemble 
after rescaling with N . 

However, this is not the case for T®^ 1 processes. Both 
the |0) and the \W) state are strongly affected by such 




FIG. 3: Collectively enhanced Rabi oscillations between an 
isolated NV qubit and a NV ensemble. The probability to 
find the qubit in the m s = 1 state, p q , goes to zero within a 
time t n ~ 600 ps. The inset shows the collectively enhanced 
coupling strength V c between the qubit and the ensemble for 
four different realizations. The dashed line shows the asymp- 
totic 1/R 3 dependence. 



processes, since the existence of any additional spin in 
the m s = 1 state corresponds to an effective magnetic 
impurity; this impurity modifies the energy of the col- 
lective state, thus tuning it out of resonance with the 
NV qubit. Additionally, this new state is also no longer 
an eigenstate of the Hamiltonian; numerical simulations 
demonstrate that this state dephases very quickly due 
to dipolar interactions within the ensemble. Thus, since 
any single spin T®^ 1 error will immediately dephase the 
collective state, the effective decoherence rate owing to 
Pj^ 1 is enhanced by N and scales with the size of the 
ensemble. 

While the system size scaling of p^ 1 errors might seem 
unfortunate, in solid-state spin systems, it is often the 
case that T\ ^> T2. Our proposed protocol is particu- 
larly useful in cases where T\ /N remains longer than T2 , 
implying that the ensemble's decoherence is dominated 
by dephasing as opposed to the enhanced depolarization. 
The specific example of NV centers highlights this crucial 
point. The dephasing of the NV originates from fluc- 
tuating mag netic fields as neighboring pairs of dipoles 
flip-flop 28, 29[. Even at low temperatures it is impos- 
sible to freeze out such magnetic fluctuations and T2 re- 
mains on the order of milliseconds 
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. On the other 
hand, the depolarization of the NV is thought to origi- 
nate from an Orbach spin-phonon process; such a process 
has an exponential dependence on temperature and im- 
plies that even moderate cooling can yield exceedi ngly 
long Ti times (^> Is at cryogenic temperatures) [3lM33|. 
By liquid nitrogen temperatures, the errors introduced 
by the enhanced T\ processes are already sub-percent, 
enabling us to focus on the effects of dephasing. An al- 
ternate approach to combat the enhanced depolarization 
of the collective state is to utilize conventional dynamical 
decoupling techniques (e.g., WAHUHA) [34| to suppress 
dipolar interactions within the individual ensembles. 

Collective quantum gates. — We now turn to a possi- 
ble application where isolated NV qubits are interspersed 
with high-density NV ensembles, forming a regular struc- 
ture, as depicted in Fig. H^i. The qubits are used for 
initialization, single-qubit rotations, and readout. Two- 
qubit gates between remote spin qubits are mediated by 
the ensemble between them and thus benefit from col- 
lectively enhanced interactions. The gate time t g is lim- 
ited by the SWAP time t n required to to transfer the 
information from one of the qubits to the ensemble (re- 
quired four times per gate operation ) jsEJ. The result- 
ing error (assuming T\/N ^> T2) of the gate is given by 
e = l — ex p[— (Atn/T^) 3 ] in the presence of spin echo 
decoupling [36|. For an error of e = 10 -2 , this translates 
to a required dephasing time of T| ff = 11 ms, which can 
be readily realized in isotopically pure diamond samples 
[H, [Soj or by using dynamical decoupling pulses 29|, HtJ . 
The requirements on the coherence time can be further 
relaxed by increasing the number of spins in the ensemble 
or by reducing the qubit-ensemble separation. 
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FIG. 4: Scalable architectures with collectively enhanced in- 
teractions, corresponding to a lattice spacing of R = 100 nm, 
compatible with sub-wavelength optical addressing, (a) In- 
dividually addressable NV qubits (yellow) are used for single 
qubit operations (SQR), whereas the collectively enhanced in- 
teraction with an ensemble is used to mediate two-qubit gates 
via SWAP operations, (b) NV ensembles are used as collec- 
tive qubits, where also single qubit operations are performed 
using the collective \W) state. 



An architecture featuring even better gate fidelities can 
be realized using a collective encoding scheme for the 
qubits (see Fig. [4)3). There, the logical |0) state corre- 
sponds to all nuclear spins being polarized, while the log- 
ical |1) state is a collective nuclear spin \W) state. This 
state can be prepared by applying a microwave pulse to 
map the electronic \W) state onto a nuclear spin \W) 
state [22]. The timescale for such a single qubit op- 
eration is limited to approximately 100 kHz by the hy- 
perfine splitting of the NV centers in the m s = 1 state 
(A\\ « -2.14 MHz for 14 N) In this collective qubit 

architecture, two-qubit gates between ensembles are en- 
hanced by a factor of N instead of a//V, thus leading to a 
SWAP time of t^ = 70 /is. Thus, we find that a gate error 
of e = 10 -2 requires a dephasing time of T| ff = 700 /is, 
while for e = 10 -4 , a coherence time of T| ff = 3 ms is 
needed [H,[30j. 

In summary, we have shown that collectively enhanced 
interactions can be realized between an NV qubit and a 
mesoscopic NV ensemble. Our proposed approach relies 
upon a transverse magnetic field to inhibit the localiza- 
tion of symmetric W-eigenstate. Our work enables the 
realization of collectively enhanced quantum gates with 
high fidelity and provides an important step towards the 
realization of scalable quantum information architectures 
involving solid-state electronic spins. 
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